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Q [ Abstract 



We study a kinetic mean- field equation for a system of particles with different sizes, in 
which particles are allowed to coagulate only if their sizes sum up to a prescribed time- 
. dependent value. We prove well-posedness of this model, study the existence of self-similar 

solutions, and analyze the large-time behavior mostly by numerical simulations. Depending 
on the parameter fcg , which controls the probability of coagulation, we observe two different 
' scenarios: For kg > 2 there exist two self-similar solutions to the mean field equation, of 

. which one is unstable. In numerical simulations we observe that for all initial data the 

rescaled solutions converge to the stable self-similar solution. For ko < 2, however, no self- 
similar behavior occurs as the solutions converge in the original variables to a limit that 
depends strongly on the initial data. We prove rigorously a corresponding statement for 
fco G (0, 1/3). Simulations for the cross-over case ko — 2 are not completely conclusive, but 
indicate that, depending on the initial data, part of the mass evolves in a self-similar fashion 
, whereas another part of the mass remains in the small particles. 
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O ■ 1 Introduction 



Mass aggregation is a fundamental process that appears in a large range of applications, such as 
formation of aerosols, polymerization, clustering of stars or ballistic aggregation, see for instance 
I [2l[6l[T2]. In all these applications, certain types of 'particles' form clusters that are characterized 

by their 'size' or 'mass' x. Smoluchowski's [Hj classical mesoscopic mean- field description of 
irreversible aggregation processes describes the evolution of the number density g{t, x) of clusters 
of size X per unit volume at time t. Clusters of size x and y can coalesce by binary collisions to 
form clusters of size x + y at a. rate given by a kernel K{x, y), such that the dynamics of g is 
given by 



oo 



d IP 

^5'(*,a;) = -y K{y,x-y)g{t,x-y)g{t,y)dy - g{t,x) J K{x,y)g{t,y) dy . (1) 



An issue of fundamental interest in the mathematical analysis of coagulation processes is the 
phenomenon of dynamic scaling for homogeneous kernels. This means that for initial data from a 
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certain class the solution to ([T]) converges to a certain self-similar solution. Unfortunately, except 
for some special kernels such as the constant one, this question is still only poorly understood 
(see e.g. [8] for an overview). While it has been common belief in the applied literature that 
self-similar solutions are unique, it has recently been shown for some special kernels [9] that there 
is a whole one-parameter family of self-similar solutions. These solutions can be distinguished 
by their tail behavior, and their respective domains of attraction are characterized by the tail 
behavior of the initial data. 

In contrast to this, very little is known for other homogeneous kernels. The existence of 
fast-decaying self-similar solutions for a range of kernels is established in [3l H], but both the 
existence of solutions with algebraic tail and the uniqueness of solutions are still open problems. 
In general, unless explicit methods such as the Laplace transform work, just proving existence 
of self-similar solutions to a coagulation equation can be a formidable task. 

Thus, despite their fundamental role, many properties of mean-field models for coagulation 
processes, in particular with respect to dynamic scaling, are not well understood. Motivated 
by an application of elasto-capillary coalescence of wetted polyester lamellas [1], we investigate 
this question for a special singular coagulation kernel K. This kernel allows only those clusters 
to coagulate that can form clusters of a given maximal size. To our knowledge this is the first 
mathematical analysis for such type of kernels. Even though tails of the size distribution do not 
play a role here, we find that self-similar solutions are still not unique and the analysis of the 
long-time behavior turns out to be delicate. 

In model considered here, two particles can coagulate only if the sum of their sizes is equal to 
M{t), where M is a given increasing function of time. This means that at time t only particles 
of size M[t) are emerging and the amount of particles of all smaller sizes is decreasing. This 
corresponds to K{x,y,t) ~ 6o{x + y — M{t)), where 6o denotes the delta-distribution in and 
the factor of proportionality may depend on time. 

As above, we denote the number density of particles of size x at time t by g{t, x). The total 
number of particles per unit volume at time t is then N{t) = jj^^*^ g{t, x)dx. At time t, the 
density of particles of any size x < M{t) is decreasing at a rate proportional to the probability 
density of a particle of size x meeting a particle of size M{t) — x. The coagulation process can 
hence be described by 



where K{t) is a rate function proportional to the expected number of coagulation events per 
unit time. Motivated by [Ij we make the ansatz K{t) = koM'(t)N(t), where fco is a constant 
of proportionality that depends on the particular physical process. The above equation hence 
reads 



The coagulation process described above neither creates nor destroys mass, which is expressed 
by the mass conservation equation 



where cr is a constant. From ([2]) and ([3]) we will derive an equivalent condition for g{t, M{t)), 
see Equation ([7]) below. 

For X > M{t), g{t,x) is not changing since particles of size greater than M{t) can not form 
a particle of size M{t) by coagulation. Hence x) = for x > M{t). We are normally 





(2) 



■Mit) 



xg{t, x)dx = a 



(3) 
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interested in processes where the value of M at the starting time is greater than the size of ah 
particles existing initially, and hence we also assume g{t,x) = for x > M{t). 

An important feature of Equation ([2]) is its invariance under reparametrization of time. As a 
consequence M{t) is not determined by the initial data but can be prescribed to be an arbitrary 
increasing function of time (see also [71 [10], where such an invariance has been crucial in the 
analysis of a model for min-driven clustering). In the following we choose M(t) = t and also 
normalize the mass by setting cr = 1. Consequently, in what follows we study the system of 
equations 



with t > 1 and < x < t. 

The first aim of this article is to establish well-posedness of the initial value problem and 
to study the long time behavior of solutions to (dj). To this end, we prove global existence and 
uniqueness of mild solutions in Section [21 Theorem [T] by rewriting ((4|) as a fixed point equation. 
In contrast to coagulation equations with more regular kernels, for which well-posedness can 
often be proved in the space of probability measures [5], here we need to work with functions 
that are continuous up to the points x = 1 and x = t. The reason that the fixed-point argument 
is not quite straightforward is that at any time there is an influx of particles of the largest size 
M(t) that leads to a nontrivial boundary condition. 

In Section [3] we study the existence of self-similar solutions. Combining rigorous arguments 
with numerical simulations, we identify ko = 2 as a critical value: For A;o < 2 no self-similar 
solutions exist but for k^ > 2 we find two self-similar solutions which have different shape and 
stability properties. Section [4] is devoted to numerical simulations of initial value problems. Our 
results for kg > 2 suggest that, after a suitable rescaling, each solution converges to the second 
self-similar solution. For kQ < 2, however, g converges to a steady state goo{x), x G M"*", whose 
shape depends strongly on the initial data. We finally prove this assertion in Proposition [6] for 



2 Global existence and uniqueness of solutions 

In this section we derive a notion of mild solutions to @ that relies on an appropritate refor- 
mulation of the mass constraint xg{t, x)dx = 1. We then use Banachs's Contraction Mapping 
Theorem to prove the local existence and uniqueness of mild solutions, and finally employ some 
a priori estimates to show that these mild solutions exist globally in time. 

2.1 Notion of mild solutions and main result 

To reformulate the mass constraint ([Us we suppose that a piecewise smooth solution x) to 
([Hi and ([11)2 is given. We then have 












and due to ([ll)i we obtain 




(6) 
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Moreover, substituting x x — t we find 



xg{t,x)x{t,t — x)dx = t I g{t,t — x)g{t,x)dx — j xg{t,t — x)g{t,x)dx, 
so ([6D implies 



g{t, x)g{t, t — x)dx. 



(7) 



This condition is equivalent to provided that g{t, x) solves and for the ease of notation 
we introduce the operator B : g ^ B[g\ by 



Bm) :-- 



ko 



g{t,x)g{t,t - x)dx, 



where depends on g via i^2- 

In what follows we fix some T > and introduce X as the space of all functions g that 
satisfy 

i) g:n = [l,T]x[0,T]^R+, 

ii) gil,x) = gini{x) for all x G [0, 1], 

iii) g{t, x) = if X > t, 

iv) g is continuous in all points {t, x) with x 1 and x ^ t, 
where the initial data gi^i are supposed to satisfy 



ffini G C{[0, 1]; / xgini{x)dx = 1. 

Jo 





To introduce the notion of mild solutions we now define the operator T : X ^ X hy 

ft 



rb](t,x) 



5fini(x)exp -k 



Co 



g[s,s- X] 



S[5r](x)exp ( -fco 




X{s) 
9{s,s- x] 
Nis) 



ds 



ds 



for X < 1, 

for 1 < X < t, 
for t < X, 



(8) 



with N{t) = Jq x)dx. Notice that T maps X into itself, and that for each g £ X the function 
g = T[g] satisfies in almost all points {t,x) the linear differential equation 



kg 

N{t) 



g{t,x)g{t,t - x) x<t, 



(9) 







X > t 



with N{t) = Jq g{t,x)dx and initial data g{l,x) = gmi{x) ■ Our main result can be summarized 
as follows. 

Theorem 1. For given initial data gi^i with ([5]) and any T > there exists a unique mild 
solution u € X to Q that satisfies ([8|) with T[g] = g. 
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2.2 Existence proof for fixed points of T 

To employ the Contraction Mapping Principle we identify a subset S C X that is invariant 
under T and a metric such that X is complete, S is closed, and F is a contraction. In what 
follows S is given by 



S=\~geX 



V(t, x)en: < g{t, x) < f{x), \ < Vdit, y)dy < ^ | , 



where / : [0, T] — t- R^" will be identified below. This set S is invariant under F provided that 
each g £ S satisfies 



o<T[g]{t,x)<f{x) y{t,x)en, 

\< I xT[g\{t,x)dx <l Vt G [l,r]. 



(10) 
(11) 



Towards (flOll we estimate 



ft 1 /■* 

N{t)= I g{t,x)dx>— / xg{t,x)dx> 
Jo T Jq 



2r' 



g{sj_s-_x) 

1 iV(s) 



c^s < 5ini(x) < /lo := llSil 



gjs^s- x) 

Nis) 



ds < Bmx) < Tko / f{y)fix - y)dy 



and this implies 1/N{t) < 2T. Moreover, 

i) if X < 1, then 

^g]{t,x) = 5(ini(x)exp ^-fcc 

ii) if t < X, then T[g]{t,x) = < /(x), 

iii) if 1 < X < t, then 

'^[g]{t,x) = B[g]{x)exp (^-k 

Therefore, / satisfies ([TO]) provided that 

fix) 

f{y)f{x - y)dy < for ah x > 1, and fix) > ho for all x < 1. (12) 

Tko 

Lemma 2. There exists a constant D > such that /(x) = hoe^^^ satisfies U^] . 
Proof. There is nothing to show for x < 1. For x > 1, condition ()12p can be rewritten as 

^> r eD{y^+{-yf-')dy = r e-^D-y(-y)dy. 
Tko Jo Jo 

Now suppose that A € [0, 1/2] is given. Then for all x > 1 we have x>x — ^>-^>0, and 
this implies 



-3Dxy{x+y) 



dy 



+ 







+ 



-3Dxy{x+y) 



dy 



X 



2A 
x 
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For the first inequality we liave estimated the integrand in the first and third part of the integral 
by 1, and in the second part by its value on the boundaries, which can be done since the integrand 
is a convex function. The second inequality then follows as A, D, and x are non-negative. 

We now choose A with A < (4:Tko)~^, and this implies ^ < 2Tko x > 1. Next 



we choose D large enough so that e' 
is decreasing for x > 1 we find xe 



-3DA 



< 



-3DAx 



< 



2Tko 
1 

2Tko 



and 3D A > 1. Since the function xe "^^^^ 
for all X > 1. Hence our choice of A and D 

□ 



guarantees (jl2]) . 
We now define 

Mt := max{/(t) : t G [0,r]}. 
Notice that g € S implies g{t,x) < Mt for all € il. 

Lemma 3. We can choose T > such that is satisfied for all c € S*. 

Proof. For g & S we have 



d_ 
di 



T[g]{t,x)git,t 



because r[^] is piecewise continuously differentiable in t and satisfies ([9]). The continuity prop- 
erties of r[^] imply that xT[g]{t,x)dx is differentiable in time, and we estimate 



jJ^xV[g\{t,x)dx = tT[g]{t,t)- ^r[~g]{t,x)dx 



ko 



< 



N{t) Jo 
ko 



xg(t, x)g{t, t — x)dx 



ko 
N{t) Jo 



xT[g]{t, x)g{t, t — x)dx 



xg{t,t- x) \g - T[g] \ {t,x)dx 



N{t) Jo 
< AkoT^M^. 

For t = 1 we have xr[^](l, x)dx = xgirii{x)dx = 1, so using the above bound we can choose 
T such that /q r[^](t, x)(ix stays between 1/2 and 3/2 for all t G [1, T]. □ 

We have now shown that S is invariant under F, and that there exists a constant Mt such 
that c{t, x) < Mt holds for all g e S and aU {t,x) € 17. 

In the next step we construct a norm for X such that F is a contraction on S. To this end, 
we define 



\9\\i ■-- 



sup \g{t,x)\, 



\9\\2 



sup 

te[i,T] Jo 



\g{t,x)\dx, 



and derive an estimate for |lF[^i] — F[^2]|li,2 in terms of ||^i — 52||i,2- Afterwards we show that 
F is a contraction with respect to some linear combination of these norms. 

Lemma 4. There exists a constant L > such that for any gi,g2 ^ S 
l|r[5i] -rfellli < L{T - l)\\gi - C2\\i + L\\gi - hh, 

||r[5l]-F[52]||2 < L(r-l)||5i-C2||l+L(T-l)||5i-g2||2. 
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Proof. Let (t, x) G Q. Since |e° — e''! < |a — 6| for all a,b > 0, we obtain the following inequalities. 
For X < 1 we find that 



i-r\~ lu \ r\~ w ^ h i. i^9i{s,s-x) g2{s,s-x) 
\T[gi]{t,x) - T[g2\{t,x)\ < /iq/co ' 



ds 



(13) 



Ni{s) N2{s) 

< hoko f I 2T|ffi s-x)-Us,s-x)\ + Mr '^l.|"\~-.^'^"^' 1 ds 

Ji [ Ni{s)N2{s) J 

< (2/ioTA;o||5i - 52II1 + 4T^MThoko\\gi - galb) ^ ds 

< 2hoTko{T - l)||5i - 52 111 + AT^MrhokoiT - l)\\gi - 52 lb- 
Now let X > 1, and set E[g](t,x) := exp(— /cq /* ^ds). Then, 

iV [s) 

\r[~gi]{t,x)-r[~g2]{t,x)\ = \B[~g^]{x)E[g^]{t,x) - B[g2]{x)E[g2]{t,x)\ 

< \B[gi]{x)-B[g2]{x)\-\E[~gi]{t,x)\ + 
\B[g2]{x)\-\E[g,]it,x)-E[g2]it,x)\ 

< \B[gi] (x) - B[g2] (x) \ + T^koM^\E[gi]{t, x) - E[g2] (t, x) \ . 
We treat the last two summands separately. Analogously to (|13p we estimate 

\B[gi]{x) - B[g2]ix)\ < (ITMrko + 2koM^T^)\\h - hh- 
On the other hand, we can estimate \E[gi]{t,x) — E[g2]{t,x)\ by 

/ {2Tko\gi{s, s-x)- 52 (s, s - x)\ + AMTT^ko\m{s) - N2{s)\}ds 

J X 

< 2Tko{T - l)\\gi - 52II1 + 4.T^ko\\gi - hh- 

Combining these results we find a constant L' that depends polynomially on kQ,T, Mt such that 

\T[gi]{t, x) - T[g2]{t,x)\ < L'{T - l)\\~g^ - 52II1 + L'{T - l)\\~gi - hh, x < 1, 
\r[g,]{t,x)-r[h]it,x)\<L\T-l)\\gi-g2\\i+L'\\gi-g2h, x > 1. 

To derive the bounds for the second norm we split the interval [0, T] = [0, 1] U [1, t] U [t, T], and 
using (|14p we find 

/ \T[gi]{t,x) - r[g2]it,x)\dx < L'T{T - l)\\~gi - gsHi + 2L' {T - l)\\gi - ^2- (15) 
JO 

With L := max{TL', 2L'} we then derive from ([HD, ^ that 

\r[gi]{t,x)-T[g2]{t,x)\ < L{T - l)\\gi - hWi + L\\9i - hh, 

[ \T[gi]{t,x) - r[g2]it,x)\dx < L{T - l)\\gi - hWi + L{T - l)\\gi - 52II2. 
Jo 

The assertions now follow by taking the supremum in the above inequalities. □ 



7 



Now let /3 > 2L, where L is as in the proof of Lemma HI and define a norm on X by 

For any gi, §2 ^ S we then have 

\\r[gi] - rfelll < L{1 + /3)(r - l)||5i - 52II1 + L{1 + (3{T - l))||5i - hh, 

and it is possible to choose T such that /3(T — 1) < 1/2 and L(l + /3)(T — 1) < 1/2. Hence 
L(l + /3(r - 1)) < 3/3/4, and this gives 

13 3 

l|r[5i] -r[52]!l < ^llsi -52II1 + ^/3||5i -52II2 < -^Wgi -hW- 

X equipped with |[-|[ is a Banach Space and S is a closed and bounded subset, so the Banach 
Fixed Point Theorem guarantees that F has a unique fixed point 5 € 5. By construction, this 
fixed point solves the differential equation for t < T. Moreover, since g{t,t) = B[g](t), it 
also satisfies condition ([7]), which is equivalent to ^3. 



It remains to prove that there exists a solution for all 1 < T < cx). This can be done by standard 
methods because (i) for each T > 1 there exists a constant D = D{T) > as in Lemma [21 and 
(ii) each solution satisfies Jg a;F[^](f, x)dx = 1 for all 1 < t < T. 



3 Self-similar solutions 

In this section we describe self-similar solutions to (jl]). These satisfy 



a(t)G (-) X <t, 
X > t, 

where G : [0, 1] M.q is a continuously differ entiable function, and the mass constraint ([Us 
requires 

ft ft />! 



1 



xg{t,x)dx= / xa{t)G (^^^ dx = t^a^t) / yG{y)dy 



This means t^a(t) is a positive constant. By rescaling a and G we can ensure that this constant 
is 1, so each mass preserving self-similar solution takes the form g{t,x) = t~'^G{x/t) with 
Jq yG{y)dy = 1. Using this relation in ([Ij), and substituting y = x/t, we get 



2G{y) + yG'{y) = ^G(y)G(l -y), N = G{y)dy. 



(16) 



We first consider a simplified problem 

2Giy) + yG'iy) = DG{y)G{l - y), (17) 

where Z? > is an arbitrary constant and we do not impose the mass constraint on G. At first, 
we notice that each solution G to (I17p provides a solution G to (I16p via 



ko = D [' Gix)dx, G{y) = . (18) 

Jo Jq xG[x)dx 

We also observe that (|17p is invariant under the scaling G ~^ AG, D ~^ X~^D. Consequently, in 
order to characterize the solution set of (|16p . we have to investigate (|17p for only one value of 
D, and then consider how the corresponding solutions transform under ()18p . 
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Lemma 5. For any given > and G(l/2) > there exists a unique positive solution of ([J 
on (0,1). 

Proof. We multiply (fT7|) by y and substitute F{y) = G{y)y'^ to obtain 

^lyi y(l-y)2 • ^'""^ 

Next, we decompose F into its odd and even parts with respect to 1/2 by setting 

F{y) = Fe{y) + Fo{y), 
F(l-y) = Fe(y)-F„(y), 

so (1191) transforms into 



D(Fe(y)^ - Fo{yf) ( 1 

^ Zj(Fe(|/)^ - FM') ^ ^ 

2y2(l_y)2 • 

Since (f20l) is locally Lipschitz for y G (0, 1), the local existence and uniqueness of solutions to the 
initial value problems is granted. In particular, for given values Fo(l/2) = and -Fe(l/2) > 
we find the smallest a £ [0, 1/2) such that there exists a unique solution to (f20]l on (a, 1 — a). 
The function F{y) = Fe{y) + Fo{y) then solves ([T9|) on (a, 1 — a), and satisfies 

Fte) = F(l/2)exp(/"^^^<i=), 

In particular, F is positive and by (|19p it is also increasing on (a, 1 — a). If a = then we are 
done. Otherwise we know that DF[1 — y)y~^{l — y)~^ is bounded for y G [1/2, 1 — a), which 
means that F does not blow-up at 1 — a. Thus, F^, Fq are also bounded on (a, 1 — a) and we 
can extend the solution to [a, 1 — a]. Due to the local existence result, we can then extend the 
solution to an interval (6, 1 — 6) with b G [0, a), which contradicts the minimality of a. Hence we 
have a = 0, and the proof is complete. □ 



subcritical solutions supercritical solutions moment of solution 




0. .5 1.0. .5 1. 0. 2. 8. 



Figure 1: Solutions G{y) to ^ with Z? = 1 for y G [0, 1] for different values of G(l/2) < 2 
(left) and 0(1/2) > 2 (center); moment N{G) = G{y)dy in dependence on G(l/2) (right). 

We next discuss some numerical ODE simulations that illustrate how the solutions of (|17|) 
with D = 1 depend on the value of G(l/2). For G(l/2) = 2 we get the trivial solution G = 2. 
From now on, we refer to solutions with G(l/2) > 2 as supercritical and to those with G(l/2) < 2 
as subcritical. These two types behave rather differently, see Figure [H which shows the solutions 
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for G{l/2) e {0.2, 0.6, 1.0, 1.4, 1.8} and G(l/2) G {2.0, 2.4, 2.8, 3.2, 3.6, 4.0}. Our numerical 
results indicate that each supercritical solution has precisely one local maximum between 1/2 
and 1 but no local minimum. On the other hand, a subcritical solution has two local maxima 
close to and 1, and a local minimum between 1/2 and 1, compare Figure O where 'variation 
near 1' refers to G{y) — G(l). 



subcritical solutions near subcritical solutions near 1 subcritical variations near 1 




0.00 .001 0.95 1.00 0.999 1.00 



Figure 2: Subcritical solutions from Figured] for y ~ and y ~ 1. 

Figure [U also shows how G{y)dy depends on G(l/2). Based on this, we conjecture that 
for every /cq > 2 there exist two solutions to p6|) . one having a subcritical and the other having 
a supercritical shape. For ko = 2 there is a trivial self-similar solution. For ko < 2 there seem to 
be no self-similar solution. To support this conjecture, we now prove that there is no self-similar 
solution for ko < 1: Integrating ([T6|) we find lim^_^i G{z) = N, so each solution to (fT6]) satisfies 
G{y) ~ as y ^ 0, which contradicts Jq G{y)dy < oo. For 1 < feg < 2, however, this 

argument does not apply but our numerical results indicate there is still no singular solution. 



4 Long-time behavior 

To investigate the long time behaviour of solutions to ([4]) by numerical simulations we derive a 
discrete "box model" which follows naturally from the physical interpretation. We assume the 
number of initial boxes Mb is sufficiently large (we choose Mb = 200 for all simulations) and 
consider the discrete times t = ej with j G N and e = 1/Mb. Moreover, we denote by G{j,i) 
the number of particles with size x G (ei — e, ei) at time t = 1 + ej, that means 



f-ei 

G{j, i) = e~^ / g(\ + ej, x)dx 

JeU-l) 



e(i-l) 

for all integers j > and i = 1, . . . , Mb + j. The disrete analogue to the evolution equation (jl]) 
is then given by 

G{j + l,i)-G{j,i) _{ - -^G{j,i)G{j,M^+j + 0<i<j + M^, 

^ 1 i > j + Mb + 1, 

where N{j) = e "^{=1^^ G{j, i) and the value of j+Mb + 1) is determined by the conservation 
of mass, i.e., 

G'(j,j + Mb + l) = | j^G{j,i)G{j,M,,+j + l-i). 
We now present our numerical results for three different values of /cq and three different sets 
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of initial data. More precisely, we consider kQ = 1,2,3 and assume that the initial data have 
Gaussian distributions with dispersion 0.3 and center at either 0.25, 0.5, or 0.75. 

For ko > 2 we expect convergence to one of the self-similar solutions. Numerical simulations 
suggest that the solution converges as t — )• oo to the supercritical self-similar solution that 
corresponds to k^. The same happens if the initial data is very close to the subcritical self- 
similar solution, and thus we can conclude that subcritical solutions are unstable. We failed to 
find a rigorous proof for this assertion, but the numerical evidence is strong. In Figure [3l we 
plot the scaled distributions after 0, 200, 1000 and 25000 steps (0, 1000, 25000 and 150000 if 
the center is at 0.25) for the three sets of initial data described above. The corresponding times 
are given by 0, 1, 5, and 125 (0, 5, 125, and 750). Along with these smooth curves we dotted 
the graph of the self-similar solution for fco = 3. As we see, the convergence is slowest when the 
center of the initial data is at 0.25. 



Center=0.25, kQ=3 Center=0.5, ko=3 Center=0.75, ko=3 




0. .5 1. 0. .5 1. 0. .5 1. 



Figure 3: Convergence to self-similar form for /cq = 3 

For ko < 2 we observe a completely different behaviour. We do not see any convergence 
in the self-similar scaling. On the contrary, the solutions apparently converge pointwise in the 
original variables to a limit that depends on the initial data. We will prove the corresponding 
statement rigorously in Proposition [6] for kQ < 1/3. Figure U] shows the unsealed distributions 
for fco = 1, the initial data is the same as before. The five smooth curves represents numerical 
distributions after 0, 200, 1000, 5000 and 25000 steps. 



Cenler=0.25, J(o=1 Center=0.5, /(q=1 Center=0.75, ka=^ 




0. 1. 2. 3. 0. 1. 2. 3. 0. 1. 2. 3. 4. 



Figure 4: Convergence in the original variables for ko = 1 

For A;o = 2 we expect convergence to the trivial self-similar solution but the numerical 
simulations do not support this assertion. We can conclude that either the long time behaviour 
is more complicated or the convergence is extremely slow. Figure [5] contains the distributions 
after 0, 1000, 5000, 25000 and 150000 steps and the dotted plot of the self-similar solution. 
In the first picture we observe a behavior similar to that for ko < 2, that means the mass is 
cumulated near the origin. It will possibly eventually disappear but this was not the case after 
any number of steps that we were able to simulate. This happens also for /cq > 2 if the initial 
data is more cummulated. For example, if we choose dispersion equal to 0.2 then 150000 steps 
is not enough to see the convergence for ko = 3 and center at 0.25. This observation is not 
surprising because if the initial data is supported on [0, 1/2), then the evolution cannot start at 
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Center=0.25, ko=2 



Center=0.5, ka=2 



Center=0.75, ko=2 




Figure 5: Evolution in self-similar variables for d = 2 



all. 

We conclude this section with a rigorous proof of a theorem supporting the conjecture that 
for ko < 2 there exists a limit function 



lim g{t,x) =: goo{x) > 0. 



(21) 



Notice that ()2ip implies that the rescaled function t'^g{t, yt) : [0, 1] — > M converges, in the sense 
of probability measures, to a Dirac measure with positive mass. 

Proposition 6. i) If there exists cq > such that N{t) > cq for all t > 1, then the relation 
(2l\) is satisfied. 

a) Ifko < 1/3 then N{t) > N{l)/2 for all t > 1. 
Proof, i) Suppose that N(t) > cq. From ([4]) and ([8]) with g = g we infer that 

g{t,x) = g{max{l,x),x) ex.p \ —ko / l^, , — ^-ds\. 



x{l,x) N{s) 



By assumption, we also have 



* g(s, s - x) ^ 1 /"* / N , 

— — ds < — / g{s, s — x)ds < oo, 

max{l,x) ^^\^) ^0 Jmax{l,x) 

SO g{t,x) is bounded from below. Since it is also decreasing in t for all x, there exists goo{x) as 
in din). 

ii) Due to dH, we have —2N{t)N'{t) = ko g{t,x)g{t,t — x)dx and by integration we obtain 
N{lf - N{tf = ko J g{s,x)g{s,s - x)dxds 

= ko / g{s,x)g{s, s — x)dsdx 

Jo Jmax(l,x) 

< ko g{max{l,x),x) / g{s,s — x)dsdx 

Jo J max(l,x) 

< ko g{max{l,x),x) / g{max{l, s), s)dsdx 

Jo Jo 

<ko{2N{l)-N{t)f, 
where the last inequality holds since 



(22) 



g{max{l, s), s)ds 



g{l, s)ds + I' g{s, s)ds = N{1) + (iV(l) 



Nit)). 



Now suppose that there exists t > 1 such that N{t) = N{l)/2. By ([22]) we get ko > 1/3 and 
using the continuity of N, we conclude that ko < 1/3 implies N{t) > N{l)/2 for all t. □ 
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